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Numerical  Modeling  of  Fire  Suppression  using  Water  Mist. 
2.  An  Optimization  Study  on  Jet  Diffusion  Flames 


1.  INTRODUCTION 

Water  as  a  means  of  fire  suppression  has  been  in  use  from  ancient  times.  The  phase-out 
of  halons  and  the  search  for  alternative  technologies  that  preserve  all  of  the  benefits  of  a  clean  total 
flooding  agent  without  the  adverse  environmental  impact  has  sparked  renewed  interest  in  water- 
mist  technology.  Fine  water-mist  relies  on  relatively  small  (less  than  200  pro)  “Class  1”  droplet 
sprays  to  extinguish  fires.  Although  the  usefulness  of  water-mist  fire  suppression  systems  has  been 
demonstrated  in  a  wide  range  of  applications  and  by  numerous  experimental  programs  (References 
[1]'  [9])i  a  widely  accepted  critical  concentration  of  water  droplets  required  to  extinguish  a  fire  is 
yet  to  be  determined.  Factors  that  contribute  to  the  success  or  failure  of  a  water-mist  system  for 
a  particular  application  include  droplet  size,  velocity,  spray  pattern  geometry,  nozzle  orientation, 
momentum  and  mixing  characteristics  of  the  spray  jet,  geometry  and  other  characteristics  of  the 
protected  area.  At  this  time  the  effect  of  these  factors  on  system  effectiveness  is  not  well  known. 
There  is  no  current  theoretical  basis  for  the  selection  of  spray  characteristics  and  other  important 
water-mist  system  parameters. 

Extinguishment  of  jet  diffusion  flames  by  the  introduction  of  water  sprays  at  the  base  of  the 
flame  has  been  studied  by  a  number  of  investigators  [1]-  [4],  These  flames  can  be  extinguished 
very  efficiently  by  this  technique  since  the  normal  air  entrainment  process  assures  that  all  the 
water  spray  added  adjacent  to  the  flame  is  actually  transported  into  the  flame.  The  water  required 
for  extinguishment  is  reported  in  terms  of  the  ratio  of  the  water  supply  rate  to  the  fuel  supply 
rate  at  extinction.  These  ratios  range  from  1.5  to  10.  While  no  systematic  evaluation  of  droplet 
sizes  on  the  water/fuel  ratio  at  extinction  has  been  performed,  the  available  data  indicate  that  the 
ratio  is  reduced  with  decreasing  droplet  size  for  laboratory  flames.  Large  scale  flames  could  be 
extinguished  with  a  water/fuel  ratio  of  1.6.  However,  ratios  of  up  to  10  may  be  required  depending 
on  the  orientation  and  geometry  of  the  spray  nozzles. 

The  available  research  in  water  spray  fire  extinguishment  was  reviewed  by  McCaffrey  [2],  Rash- 
bash  [5]  and  by  Tatem  et  al.  [6].  McCaffrey  [2]  has  reviewed  the  available  data  for  the  application 
of  water  sprays  as  a  jet  diffusion  flame  suppression/extinguishment  agent.  Reduction  in  flame  tem¬ 
perature  during  sub-extinguishing  application  rates  of  water  spray  has  been  observed  to  correlate 
with  a  single  spray  parameter  -  the  median  drop  diameter.  The  effect  of  adding  water  spray  to 
the  gas  flow  below  the  base  of  a  lifted  flame  is  to  change  the  flame  shape  and  to  lower  peak  flame 
temperature  and  radiation  levels.  Extinguishment  near  blow-off  was  observed  to  be  due  to  a  shift 
in  flame  position  [2].  The  influence  of  droplet  diameter,  density  and  velocity  on  flame  suppression 
for  various  nozzle  orientations  and  spray  pattern  geometries  has  not  been  reported  in  a  detailed 
parametric  fashion  in  the  literature. 
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There  has  been  some  theoretical  and  experimental  work  ( [10]-  [15])  to  develop  an  understanding 
of  gas  phase  extinguishment  of  fires.  Numerical  studies  for  a  counterflow  methane/air  diffusion  flame 
with  heptane  droplets  added  to  the  fuel  stream  have  been  reported  by  Chen  et  al.  [14].  Atreya 
[15]  has  observed  enhancement  of  the  burning  rate  at  low  water  addition  rates  followed  by  flame 
extinguishment  at  higher  application  rates  for  counterflow  diffusion  flames.  Mawhinney  et  al.  [3] 
has  also  reported  evidence  of  invigoration  of  combustion  by  the  introduction  of  water-mist.  The 
time  to  extinguish  a  fire  was  reported  as  a  function  of  fuel  type,  geometry  and  mist  characteristics. 
Hoffman  and  Galea  [16],  [17]  have  used  the  field  modeling  technique  to  include  two-phase  fire- 
sprinkler  scenarios.  In  these  studies  the  fire  is  seen  simply  as  a  heat  source  without  the  added 
complication  of  combustion.  The  model  results  are  compared  with  experimental  data  and  shows 
good  agreement  near  the  sprinkler  source,  but  deteriorates  in  the  far  field.  More  recently,  Prasad  et 
al.  [18],  [19]  have  studied  the  relative  contribution  of  various  flame  suppression  mechanisms  (thermal 
cooling  and  oxygen  dilution),  and  have  made  detailed  comparisons  with  experimental  results. 

The  overall  objective  of  this  study  is  to  develop  and  apply  a  numerical  model  for  studying  the 
underlying  processes  involved  in  suppression  of  fires  with  water-mist.  Overall  flame  suppression 
and  mist  entrainment  into  the  flame  and  their  dependence  on  droplet  diameter,  spray  velocity  and 
injection  characteristics  (orientation  of  the  nozzles  relative  to  the  flame  and  angle  of  injection) 
is  determined.  Sectional  trajectories  are  used  to  identify  the  regions  of  the  flame  (preheat  zone, 
combustion  zone  or  plume  zone)  where  the  droplets  evaporate  and  absorb  energy.  Numerical 
simulations  are  performed  for  both  symmetric  and  asymmetric  spray  injection  geometries  (base 
injection  or  side  injection).  In  this  report,  “Flame  Suppression”  is  used  to  indicate  a  reduction  in 
the  net  heat  output  from  a  flame.  The  water-mist  required  for  different  levels  of  suppression  is 
reported  in  terms  of  the  ratio  of  water-mist  flow  rate  to  the  fuel  flow  rate.  This  ratio  is  related  to 
the  water-mist  loading  which  is  simply  the  ratio  of  water-mist  flow  rate  to  the  air  flow  rate.  The 
goal  of  the  current  work  is  not  to  optimize  the  design  of  water-mist  systems,  but  to  determine  the 
optimum  injection  characteristics  of  a  water-mist  system  to  suppress  a  specific  fire  -  a  jet  diffusion 
flame.  Several  such  fire  science  studies  are  needed  before  arriving  at  the  optimum  characteristics 
for  any  arbitrary  fire.  Additional  system  configuration  studies  will  be  needed  to  determine  the 
optimum  placement  of  water-mist  systems  for  suppression  of  fires  in  a  specific  enclosure. 

2.  PROBLEM  FORMULATION 

We  consider  the  interaction  of  fine  water-mist  with  methane-air  diffusion  flames  stabilized 
above  a  Wolfhard-Parker  burner  geometry.  There  are  two  interacting  physical  phases,  the  gas 
phase  involving  convection,  diffusion  and  chemical  reactions  between  the  constituent  species,  and 
the  liquid  phase,  representing  the  evaporating  water  droplets.  A  two-continuum  formulation  is 
adopted  in  which  gas  properties  and  the  droplet  properties  are  each  described  by  equations  of  the 
Eulerian  form. 

2.1  Gas  Phase  Equations 

Detailed  modeling  of  gas-phase  reactive  flows  is  based  on  a  generally  accepted  set  of  time- 
dependent  coupled  partial  differential  equations  maintaining  conservation  of  total  mass,  momentum, 
total  energy  and  individual  species  density.  A  strong  conservation  form  of  the  two-dimensional, 
unsteady,  compressible  Navier-Stokes  equations,  used  to  describe  gas  phase  reactive  flow  systems 
can  be  written  as  follows  [18] 
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In  these  equations  x  and  y  denote  the  independent  spatial  coordinate  and  t  denotes  the  temporal 
coordinate;  p  the  mass  density;  pk  the  density  of  the  kth  species;  u  and  v  are  the  x  and  y  components 
of  the  fluid  bulk  velocity;  P,  the  pressure;  E ,  the  total  energy  of  the  fluid  per  unit  mass;  g,  the 
gravitational  body  force  per  unit  mass;  Uk  and  14  are  the  two  components  of  the  diffusion  velocity 
for  the  kth  species;  txx,  rxy,  ryy  are  the  components  of  the  stress  tensor  for  newtonian  fluid  in 
rectangular  coordinate;  qx  and  qy  are  the  *  and  y  components  of  the  heat-flux  vector.  mevap  is  the 
source  term  arising  in  the  mass  conservation  equation  due  to  evaporation  of  the  water  droplets, 
resulting  in  the  production  of  water  vapor.  The  corresponding  term  in  the  water  vapor  species 
conservation  equation  is  represented  by  ihk,eVap ♦  <^k  is  the  rate  of  production  of  the  kth  species 
due  to  chemical  reactions.  Qchem.  is  an  exothermic  source  term  in  the  energy  equation  representing 
chemical  energy  release,  whereas  Qevap  represents  absorption  of  energy  due  to  droplet  evaporation; 
Qrad  is  the  radiative  heat  loss  term.  Sx  and  Sy  are  the  source  terms  in  the  x  and  y  momentum 
equations  representing  the  cumulative  drag  force  exerted  by  the  droplets  on  the  gas  phase. 


Ideally,  we  would  like  to  simulate  the  chemical  reactions  by  including  a  detailed  set  of  elemen¬ 
tary  reactions  to  describe  the  production  of  the  individual  species  and  the  energy  release  in  the 
flame.  The  cost  of  computer  time  and  memory  required  to  track  the  individual  species  makes  this 
prohibitive  for  problems  in  which  parametric  studies  are  planned.  Therefore,  the  chemical  reaction 
and  energy-release  process  for  methane-air  combustion  is  described  phenomenologically  based  on  a 
single  step  reaction, 


CH4  +  202  -4  C02  +  2H20, 


(6) 


using  a  finite- rate,  quasi-global  Arrhenius  expression.  Westbrook  et  al.  [20]  has  developed  simpli¬ 
fied  reaction  mechanisms  for  the  oxidation  of  hydrocarbons  using  a  laminar  flame  model  and  has 
prescribed  an  overall  rate  of  consumption  of  methane  as 


<j[0ff4] 

dt 


-1.3^8exp(-484OO/f?T)[0if4]-o-3[O2]1-3. 


(7) 


The  above  expression  is  used  to  compute  the  depletion  of  methane.  The  various  diffusion  coefficients 
are  obtained  from  detailed  kinetic  theory  and  are  fitted  over  a  suitable  temperature  range  using  a 
third  order  polynomial  [21]. 
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2.2  Sectional  Water  Mist  Model 


In  this  section  the  equations  describing  a  vaporizing  spray  are  obtained  by  a  standard 
control- volume  approach.  It  is  assumed  that  coalescence  and  breakup  of  droplets  are  insignificant 
with  in  the  computational  domain.  An  Eulerian  approach  is  adopted  and  the  droplet  properties  are 
treated  as  if  they  were  continuous  in  the  domain  as  the  gas  properties  [22],  [23].  The  vaporization 
process  can  be  described  by  a  set  of  coupled  differential  equations  for  the  concentration  n,-  of  the 
discrete  droplet  sizes  per  unit  volume  of  the  fluid, 

dn{  ®ni  dfl{  _  _  ,  . 

f*  ^7,2 I-  —  -EiUi  +  (8) 

where  u\j  and  vij  are  the  x  and  y  components  of  the  velocity  vector  of  the  ith  droplet  size  and  E{ 
is  the  frequency  of  the  molecule  evaporation.  This  frequency  appearing  in  the  discrete  form  of  the 
droplet  population  balance  equations,  depends  on  surface  area  of  the  droplet,  temperature  of  the 
ambient  gas  and  other  flow  parameters.  Functional  forms  of  this  frequency  can  be  obtained  through 
known  expressions  of  vaporization  rates  of  single  droplets  or  may  be  determined  through  correlations 
with  known  experimental  results  of  measured  vaporization  rates.  Since  the  total  number  of  droplet 
sizes  needed  to  simulate  actual  fuel  sprays  can  be  immense,  sectional  conservation  equations  have 
been  employed  [24],  The  method  of  “section  conservation”  avoids  the  dimensionality  problem 
associated  with  the  discrete  form  of  the  droplet  population  balance  equations  8.  This  method  is 
based  on  dividing  the  droplet  size  domain  into  sections  and  dealing  only  with  one  integral  quantity  in 
each  section.  This  sectional  representation  has  the  advantage  that  the  integral  quantity  is  conserved 
within  the  computational  domain  and  the  number  of  conservation  equations  is  substantially  reduced 
so  as  to  equal  the  number  of  sections. 


We  divide  the  entire  droplet  size  domain  into  M  arbitrary  sections,  and  define  Qj  to  be  an 
integral  quantity  of  the  spray  within  the  jth  section.  Thus 

rvi 

Qj=  /  vndv,  j  =  1,2,  (9) 

Jvj-l 

where  n(v,t)  is  the  number  concentration  function  and  Vj-\  and  Vj  denote  the  volumes  of  the 
smallest  and  largest  droplets,  respectively,  in  section  j.  The  quantity  Qj  can  be  related  to  volume 
or  the  mass  of  the  droplets  within  the  jth  section.  Using  equations  8  and  9,  the  sectional  conservation 
equation  can  be  written  as 

+  +  =  +  i  =  l,2,...,Af,  (10) 

where  Bjj+i  and  Cj  are  the  sectional  coefficients.  Assuming  a  d?  law  to  re-express  the  sectional 
vaporization  coefficients  in  terms  of  droplet  diameters  [18]  and  assuming  a  continuous  division  into 
sections  the  coefficients  can  be  expressed  as 


BjJ+1  =  -E(T) 


Cj  =  2  E{T) 


dL, 


J+l 


ldU,j+l  ~  dL,j+ 1  J 


3  —  1, 2, . . .,  M  —  1, 


dL,j 


3{dUJ  ~  dL,j) 


L  dU,j  -  dL,j 


dU,j~dL,j  J 


,  j  =  l,2,...,Af, 


(11) 
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where  E(T )  is  the  burning  rate  coefficient,  which  depends  on  temperature  and  other  properties  of 
the  droplet  and  its  surroundings.  For  example,  with  an  initial  injection  droplet  diameter  of  150 y, 
the  droplet  size  domain  is  divided  into  5  continuous  division  (M  =  5).  These  sections  are  0  —  30^m, 
30  -  60 ym,  60  —  90 ym,  90  -  120 pm  and  120  —  150/zm. 


Each  of  the  droplet  sections  is  assumed  to  have  its  own  unique  velocity  different  from  that  of 
the  gas  phase.  Momentum  conservation  equations  are  formulated  for  each  droplet  section  and  are 
coupled  to  those  of  the  gas  phase  through  the  phase  interaction  terms  (drag  terms).  The  droplet 
drag  equations  are  given  by 


duij 

8t 


+  uij 


duU 

dx 


+  v‘,j 


duij 

dy 


dvU 

dt 


+  ui,j 


dvU 

dx 


+ 


dvij_ 

dy 


—  Fx,j  i  j 


—  Fyj  5) 


=  1,2 AT, 
j  =  1,2,  ...,M, 


(12) 


where  Fxj  and  Fyj  are  the  x  and  y  components  of  the  drag  force  acting  on  the  jth  sectional  velocity 
modeled  based  on  Stokes  Law.  The  interphase  drag  force  experienced  by  the  two  phases  is  related 
to  the  particle  Reynolds  number  by  a  drag  coefficient  Cd •  The  gas  phase  momentum  source  terms 
Sx  and  Sy  are  equal  and  opposite  to  the  sectional  density  weighted  sum  of  Fx  j  and  Fy  j  respectively 
[18],  [25]. 


2.3  Numerical  Algorithm 

The  governing  equations  are  rewritten  in  terms  of  finite-volume  approximations  on  an 
Eulerian  mesh  and  solved  numerically  using  a  time  marching  procedure.  The  fluid  convection 
is  solved  with  a  high-order  implicit  algorithm,  Barely  Implicit  Correction  to  the  Flux-Corrected 
Transport  (BIC-FCT)  [26],  that  was  developed  to  solve  the  convection  equations  for  low- velocity 
flows.  The  Flux-Corrected  Transport  (FCT)  [27]  algorithm  itself  is  an  explicit,  finite-difference 
algorithm  that  is  constructed  to  have  fourth-order  phase  accuracy.  The  solution  approach  consists  of 
separate  algorithms  for  each  of  the  individual  processes,  which  are  coupled  together  by  the  method 
of  fractional  time-step  splitting.  The  algorithms  for  convection,  thermal  conduction,  molecular 
diffusion,  chemical  reactions  and  the  coupling  of  the  individual  processes  have  been  previously 
discussed  in  detail  [26]. 

Figure  1  shows  a  schematic  diagram  of  a  2-Dimensional  Wolfhard-Parker  diffusion  flame  burner 
along  with  the  associated  computational  domain  in  which  the  solutions  are  desired.  The  various 
water-mist  injection  configurations  such  as  base  injection,  side  injection  and  top  injection  used  in 
this  optimization  study  have  also  been  illustrated.  A  specific  solution  of  the  reactive  flow  equations 
is  determined  by  the  initial  and  boundary  conditions  that  describe  the  geometry  of  the  system  and 
exchange  of  mass,  momentum  and  energy  occurring  between  the  system  and  the  rest  of  the  physical 
world.  The  various  species  densities  and  jet  velocities  are  prescribed  at  the  inflow  boundary.  At 
the  outflow  boundary,  the  normal  gradients  of  total  density,  species  densities  and  momentum  was 
assumed  to  be  zero.  The  pressure  at  large  distances  from  the  burner  surface  is  assumed  to  be  equal 
to  the  ambient  pressure.  Symmetric  boundary  condition  was  employed  at  the  center  line  of  the 
computational  domain  by  use  of  anti-symmetric  reflection  of  tangential  velocity  v,  and  symmetric 
reflection  of  all  other  variables.  A  slip  wall  boundary  condition  is  employed  at  the  lateral  boundary 
by  assuming  that  the  flux  of  all  transported  properties  across  the  wall  is  zero.  For  asymmetric 
flame  calculation  slip  wall  boundary  conditions  are  invoked  at  both  lateral  boundaries. 
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The  computational  domain  for  a  symmetric  calculation  consists  of  a  stretched  96  x  128  grid, 
concentrated  near  the  contact  point  of  the  fuel  and  the  oxidizer.  The  smallest  grid  cell  is  a  0.5mmx 
0.5mm  square  cell.  Numerical  simulations  were  performed  by  halving  the  grid  size  to  check  for 
accuracy  and  consistency.  The  algorithms  have  been  optimized  for  efficient  performance  on  vector 
machines.  Typical  computational  time  on  a  C-90  single  processor  machine  was  approximately  4 
hours  to  reach  a  steady  state  solution. 

3.  RESULTS  AND  DISCUSSION 

The  governing  equations  described  earlier  were  solved  using  the  numerical  models  briefly 
discussed  above,  to  simulate  co-flow  methane-air  diffusion  flames  stabilized  above  a  Wolfhard  Parker 
burner.  Computations  were  performed  to  study  the  interaction  of  water-mist  and  the  diffusion  flame 
for  various  droplet  diameters,  spray  injection  density,  velocity  and  injection  configurations  (side 
injection  or  base  injection,  symmetric  or  asymmetric  spray  pattern  geometry).  Simulations  were 
performed  to  estimate  the  impact  of  various  injection  configurations  on  overall  flame  suppression, 
flame  structure  and  flow  pattern. 

First,  numerical  simulations  were  performed  for  the  base  case  methane-air  diffusion  flame  (with¬ 
out  any  injection  of  water-mist).  The  average  exit  velocity  of  fuel  jet  (CH4)  was  2.81  cm/s  and  the 
average  exit  velocity  of  the  air  jet  was  18.1  cm/s  (See  Figure  1).  Temperature  profiles  at  various 
heights  above  the  burner  surface  were  compared  with  experimental  data  obtained  over  a  simi¬ 
lar  geometry.  At  each  height  the  numerically  computed  temperature  profiles  compared  favorably 
with  the  thermocouple  temperature  measurements  [18],  indicating  that  the  exothermic  chemical 
reactions  and  energy  release  are  adequately  modeled  by  the  single  step  Arrhenius  reaction  mecha¬ 
nism.  Numerical  calculations  and  temperature  measurements  were  performed  for  various  nitrogen 
dilution  levels  on  methane-air  diffusion  flames  to  study  the  impact  of  oxygen  displacement.  Nu¬ 
merically  computed  temperature  profiles  again  compared  favorably  with  experimental  results  at 
various  heights  above  the  burner  surface  [18],  [19]. 

3.1  Base  Injection  Configuration  -  Co-flow  Injection 

With  the  confidence  gained  by  the  comparison  with  experimental  data,  a  series  of  numeri¬ 
cal  simulations  were  performed  to  determine  the  underlying  processes  involved  in  suppression  of 
flames  with  water-mist.  Water-mist  is  introduced  along  with  the  air  co-flow  at  the  base  of  the 
diffusion  flame  (base  injection).  Parametric  studies  were  performed  with  50 pm  droplets  injected 
uniformly  over  the  length  of  the  air  channel  with  spray  injection  density  varying  from  1000  dr  op/ cm3 
(0.005kg/m3)  to  4000  drop/ cm?  and  injection  velocity  varying  from  25  cm/s  to  100  cm/s.  Similar 
studies  were  performed  for  100pm  and  150 pm  monodisperse  droplets  with  injection  density  vary¬ 
ing  from  50  drop/cm3  to  1000  drop/cm3  and  mist  injection  velocity  varying  over  the  range  of  25 
cm/s  to  500  cm/s.  The  temperature  and  heat  release  rate  profiles  indicate  an  increase  in  flame 
height  and  flame  suppression  due  to  the  sub-extinguishing  application  rates  of  water-mist.  The 
heat  release  rate  profile  indicates  the  presence  of  an  endothermic  heat  release  rate  profile  (negative 
value)  showing  the  approximate  location  of  evaporation  of  the  droplets.  These  profiles  along  with 
the  section  droplet  density  contours  indicate  that  droplets  are  unable  to  cross  the  flame  sheet  and 
enter  the  diffusion  flame. 

The  net  suppression  due  to  water-mist  and  the  impact  of  droplet  diameter,  injection  velocity 
and  spray  density  have  been  summarized  in  Figure  2.  The  abcissa  shows  the  net  amount  of  water- 
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mist  injected  per  unit  mass  flow  rate  of  fuel  (mmist/ihfuei),  whereas  the  ordinate  is  a  ratio  of 
the  integrated  heat  release  in  the  computational  domain  with  water-mist  (A Hj)  to  the  integrated 
value  without  water-mist  (A H0).  This  ratio  therefore  represents  a  net  integrated  measure  of  the 
suppression  of  a  flame.  Efficient  design  of  water-mist  systems  aims  at  obtaining  the  maximum 
amount  of  suppression  (values  of  AHt/AH0  close  to  0.0)  with  the  minimum  amount  of  water- 
mist  added  to  system,  that  is,  to  minimize  the  value  of  the  abcissa  and  the  ordinate.  Along  each 
curve  in  Figure  2,  droplet  diameter  and  injection  velocity  is  held  constant,  but  spray  density  is 
progressively  increased  resulting  in  higher  mist  flow  rates.  For  each  droplet  diameter,  the  net 
suppression  increases  with  spray  density.  The  ratio  of  mist  to  fuel  flow  rate  ( ihmist/ihfuel )  required 
for  extinction  of  the  diffusion  flame  is  approximately  10  for  50 pm  droplets  and  increases  to  40  for 
150 pm  droplets.  Figure  2  also  shows  cases  with  three  different  spray  injection  velocities  (25 cm/s, 
250cm/s  and  500cto/.s).  As  the  spray  velocity  increases,  we  find  that  net  suppression  reduces  for 
a  given  (riimist/mfuei)  ratio  due  to  reduced  entrainment.  Overall  our  results  indicate  that  for  the 
co-flow  configuration,  smaller  diameter  droplets  produce  maximum  suppression  for  the  minimum 
spray  mass  density. 

3.2  Symmetric/Asymmetric  Injection  Configuration 

Computations  were  carried  out  with  water-mist  injection  directed  at  an  angle  of  45°  to  the 
burner  surface  for  a  wide  range  of  droplet  diameters,  injection  densities  and  spray  velocities.  Both 
symmetric  and  asymmetric  spray  geometries  were  considered.  Symmetric  injection  refers  to  water- 
mist  injection  that  is  directed  towards  the  diffusion  flame  from  both  sides  of  the  flame  (Figure 
3).  Under  these  conditions,  the  symmetric  boundary  condition  can  be  invoked  at  the  centerline. 
Asymmetric  injection  refers  to  mist  injection  directed  towards  the  flame  on  one  side  and  away 
from  the  flame  on  the  other  side.  Under  these  conditions  the  cumulative  drag  force  exerted  by  the 
droplets  on  the  gases  may  be  large  enough  to  bend  the  flame  in  one  direction.  Symmetry  conditions 
cannot  be  invoked  at  the  geometric  center  line  of  the  system  and  the  entire  computational  domain 
has  to  be  solved. 

Figure  3  shows  the  net  suppression  effect  of  water-mist  due  to  symmetric  and  asymmetric  injec¬ 
tion  at  45°  angle  to  the  burner  surface.  Water-mist  is  injected  along  with  the  air  co-flow  at  the  base 
of  the  diffusion  flame  (base  injection  configuration).  Along  each  curve  the  droplet  injection  density 
is  kept  constant  but  spray  velocity  is  changed  resulting  in  higher  mist  flow  rates.  Results  show  that 
for  50  pm  droplets  under  symmetric  injection  condition,  higher  injection  densities  (2000  dr  op /cm3) 
produce  better  suppression  as  compared  to  lower  injection  densities  (1000  drop/ cm3).  However, 
asymmetric  injection  conditions  produce  lower  suppression  even  with  higher  mass  injection  densi¬ 
ties  (4000  drop /cm 3)  as  compared  to  the  symmetric  case.  During  asymmetric  water-mist  injection, 
smaller  amount  of  water-mist  is  entrained  into  the  diffusion  flame  as  compared  to  the  symmetric 
case  resulting  in  lower  suppression.  Results  also  indicate  that  under  asymmetric  conditions,  150 pm 
droplets  result  in  better  flame  suppression  as  compared  to  50 pm  droplets.  This  result  is  in  contrast 
with  that  shown  in  Figure  2  which  shows  that  smaller  droplets  produce  better  flame  suppression. 
In-addition,  asymmetric  flames  require  higher  mmist/mfuei  ratio  for  suppression.  The  temperature 
and  streamline  pattern  indicate  that  flame  structure  and  flow  field  pattern  is  significantly  different 
between  the  symmetric  and  asymmetric  configurations. 

Temperature  contours  during  symmetric  injection  of  water-mist  at  45°  angle  to  the  air  co-flow 
have  been  shown  in  Figure  4.  The  mist  injection  velocity  is  varied  from  0.25m/s  to  2.25 m/s.  The 
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arrows  indicate  the  initial  injection  direction  at  the  base  of  the  diffusion  flame  burner.  Results 
indicate  that  as  the  mist  velocity  increases  the  flame  height  (defined  as  the  height  above  the  burner 
surface  along  the  centerline,  where  the  temperature  reaches  a  local  maxima)  increases  but  the 
maximum  temperature  decreases.  For  water-mist  injection  velocity  of  0.25m/s,  we  observe  a  peak 
flame  temperature  of  1870A'  and  a  flame  height  of  7.0cm.  On  the  other  hand  for  a  mist  injection 
velocity  of  2.25 m/s  a  peak  flame  temperature  of  1580.fi'  was  obtained  and  the  corresponding  flame 
height  was  approximately  10cm.  The  gas  phase  velocity  vectors  (color)  and  streamline  pattern 
(black)  during  symmetric  injection  of  water-mist  at  45°  angle  to  the  air  co-flow,  for  various  mist 
injection  velocities  has  been  illustrated  in  Figure  5.  As  the  mist  injection  velocity  increases,  the 
drag  force  exerted  by  the  water-mist  on  the  air  flow  also  increases  resulting  in  the  bending  of  the 
streamlines  towards  the  diffusion  flame.  The  contraction  of  the  streamtube  is  accompanied  with 
increased  entrainment  of  the  air  into  the  diffusion  flame.  This  leads  to  taller  diffusion  flames  as 
shown  in  Figure  4. 

The  temperature  contours  during  asymmetric  injection  of  water-mist  at  45°  angle  to  the  air 
co-flow  have  been  shown  in  Figure  6.  The  water-mist  droplet  diameter  is  50/tm  and  the  injection 
density  is  4000 drops/ cm3.  The  mist  injection  velocity  is  varied  from  1.75m / s  to  7.00m/s.  The 
arrows  indicate  the  initial  injection  direction  of  water-mist  at  the  base  of  the  diffusion  flame  burner. 
Temperature  contours  clearly  indicate  that  the  diffusion  flame  bends  in  the  direction  of  the  mist 
flow.  We  also  observe  that  the  peak  temperature  increases  from  1970A'  (water-mist  injection 
velocity  1.75m/s)  to  approximately  2030 K  for  a  water-mist  injection  velocity  of  7.00m/s.  The 
large  mist  injection  velocity  can  significantly  alter  the  air  flow  patterns  resulting  in  higher  peak 
flame  temperatures.  Figure  7  shows  the  sectional  density  contours  above  a  methane-air  diffusion 
flame  burner  for  asymmetrical  injection  of  water-mist  at  45°  angle.  The  sectional  density  contours 
indicate  the  approximate  location  where  the  droplets  evaporate  and  absorb  energy  from  the  diffusion 
flame.  The  asymmetrical  sectional  density  contours  for  the  0  - 10 pm  and  20  -  30/m  droplets  imply 
larger  extraction  of  heat  on  one  side  of  the  diffusion  flame  resulting  in  asymmetric  flames.  This 
explains  the  observed  temperature  distribution  shown  in  Figure  6  with  one  side  of  the  diffusion 
flame  being  cooler  than  the  other  side. 

3.3  Base  Injection  Configuration  -  Angle  Injection 

Figure  8  shows  the  overall  flame  suppression  as  a  function  of  spray  injection  angle  for  various 
spray  injection  density  and  velocity.  Both  symmetric  and  asymmetric  conditions  are  shown.  All 
the  data  points  on  any  one  curve  in  Figure  8  have  the  same  ratio  of  the  mist  flow  rate  to  fuel 
flow  rate.  Hence  the  use  of  injection  angle  as  the  variable  on  the  abcissa  instead  of  rhmtst/rh/ue;. 
Results  indicate  that  for  certain  conditions  the  net  suppression  is  independent  of  the  injection  angle. 
When  the  droplet  injection  velocity  is  a  close  to  the  gas  inlet  velocity,  the  drag  force  exerted  by 
the  gas  on  the  droplets  is  strong  enough  to  make  the  water  droplets  quickly  follow  the  gas  phase. 
As  the  spray  injection  velocity  increases  and  becomes  significantly  different  from  the  gas  phase, 
results  indicate  an  optimum  angle  for  flame  suppression.  Results  for  asymmetric  injection  of  mist 
(50/m,  density=1000  dr  op/ cm3  (0.065 kg/m3),  velocity=1.75  m/s)  show  an  optimum  throw  angle 
of  45°.  Similar  results  were  observed  for  150/m  droplets  with  high  injection  velocities  of  1.0  m/s 
injected  in  symmetric  or  asymmetric  manner.  Again,  as  expected,  suppression  under  symmetric 
conditions  was  better  than  under  asymmetric  conditions  for  similar  injection  parameters  (150/m, 
density=100  drop/cm 3  (0.17 6kg/m3),  velocity=1.0  m/s).  Figure  9  shows  temperature  contours 
during  the  symmetric  injection  configuration  of  water-mist  at  various  injection  angles.  In  these 
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simulations  the  droplet  diameter  is  50/mi,  injection  density  is  100 drops/ cm3  and  injection  velocity 
is  l.OOra/s.  The  maximum  temperatures  were  1740/f,  1720 K  and  1830 K  for  injection  angles  of 
15°,  45°  and  75°  respectively. 


3.4  Side  Injection  Configuration 

The  calculations  described  thus  far  (Figures  2-8)  are  for  water-mist  injection  at  the  base  of  the 
diffusion  flame  (base  injection  configuration).  We  now  describe  water-mist  injection  through  the 
side  wall.  Side  injection  configuration  (see  figure  1)  are  important  because  they  will  have  a  direct 
influence  on  the  entrainment  rate  of  water-mist  into  the  flame.  When  water-mist  is  injected  only 
through  one  side  wall,  the  cumulative  drag  force  exerted  by  the  droplet  sections  on  the  gas  phase 
may  result  in  bending  of  the  diffusion  flame  due  to  the  cross  flow.  Figure  10  shows  the  heat  release 
profiles  (J/m?/s)  during  uniform  injection  of  150 [im  droplets  distributed  evenly  through  the  side 
wall.  The  red  contours  indicate  the  region  of  high  exothermic  chemical  activity  at  the  flame  sheet. 
The  blue  contours  mark  the  location  where  water  droplets  evaporate  and  absorb  energy.  Unlike 
the  case  for  base  injection  of  water-mist  through  the  burner  surface,  results  clearly  show  that  water 
droplets  are  able  to  penetrate  through  the  flame  sheet  and  absorb  energy  inside  the  flame  envelope. 
Figure  11  shows  the  gas  phase  velocity  vectors  superimposed  on  temperature  contours  for  the  case 
of  side  injection  of  water-mist  into  a  diffusion  flame.  Results  clearly  illustrate  the  bending  of  the 
temperature  contours  representing  a  diffusion  flame  stabilized  in  a  mild  cross  flow.  The  change 
in  gas  phase  flow  field  due  to  the  side  injection  of  water-mist  is  clearly  illustrated  by  this  figure. 
Sectional  density  contours  (refer  Figure  12)  also  indicate  the  presence  of  0  —  10 \im  and  20  —  30 fim 
droplet  sections  inside  the  flame  sheet.  These  contours  indicate  the  locations  where  water-mist 
evaporate  and  absorb  energy  from  the  diffusion  flame. 

Figure  13  shows  the  net  integrated  suppression  during  side  injection  of  water-mist  for  various 
injection  spray  densities  and  velocities.  The  effects  of  droplet  diameter  on  suppression  has  been 
shown  as  a  function  of  mist  to  fuel  flow  rate  for  asymmetric  flames.  Results  clearly  show  that  150 fim 
droplets  produce  more  flame  suppression  as  compared  to  50 \im  droplets  for  a  given  mmist/ihfuei 
ratio.  In  contrast,  base  injection  configuration  show  that  small  droplet  diameters  (50/zm)  suppress 
flames  more  than  larger  droplet  diameters  (150 (im).  The  small  droplets  exhibit  small  characteristic 
time  for  decrease  of  relative  velocity  between  the  gas  phase  and  water  droplets.  As  a  result,  these 
droplets  quickly  follow  the  gas  phase.  Larger  droplets  tend  to  travel  with  their  injection  velocity 
and  exhibit  large  characteristic  time  for  decrease  of  relative  velocity.  When  50 jim  drops  are  injected 
with  moderate  injection  velocity  (1  m/s )  through  the  side  wall,  they  are  unable  to  cross  the  width 
of  the  air  channel  and  reach  the  diffusion  flame.  These  droplets  are  dragged  by  the  air  flow  out 
of  the  computational  domain  before  they  reach  the  flame.  Only  very  high  injection  velocities  (10 
m/ s )  for  the  50 fim  droplets  are  able  to  reach  the  diffusion  flame.  As  a  result  the  suppression  curve 
for  these  drops  remains  flat  (no  suppression)  for  mmi$t/mfuei  ratio  of  up  to  60.  The  150 fim  drops 
are  able  to  cross  the  width  of  the  air  channel,  entrain  into  the  diffusion  flame  and  provide  flame 
suppression  for  smaller  injection  velocities. 

Figure  13  also  shows  that  for  the  side  injection  configuration  lower  spray  injection  density 
results  in  more  suppression  for  a  given  mist  flow  rate.  This  is  because  to  maintain  a  given  mist 
flow  rate,  the  lower  spray  density  case  has  a  higher  mist  velocity.  Flame  suppression  increases 
with  higher  mist  velocity  for  a  given  injection  density  and  droplet  diameter.  On  the  other  hand, 
as  discussed  earlier,  under  the  base  injection  configuration,  a  higher  spray  injection  density  results 
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in  more  flame  suppression.  Flame  suppression  decreases  with  higher  mist  velocity  for  a  given  mist 
flow  rate  and  droplet  diameter. 

3.5  Top  Injection  Configuration 

Computations  were  carried  out  with  water-mist  injection  from  the  top  of  the  diffusion  flame 
(top  injection  configuration),  as  shown  in  Figure  1.  Water-mist  was  injected  through  the  outflow 
boundary  and  directed  downwards  to  the  diffusion  flame.  In  the  present  calculations  the  outflow 
boundary  was  approximately  30cm  above  the  burner  surface.  Numerical  simulations  were  performed 
with  water-mist  injection  velocity  ranging  from  20 cm/s  to  200cm/s  and  an  initial  droplet  diameter 
ranging  from  50/m  to  150/m.  For  small  initial  injection  velocity  ranging  from  20  -  50 cm/s  the 
droplets  were  unable  to  overcome  the  drag  force  exerted  by  the  hot  plume  gases.  The  droplets  were 
observed  to  quickly  reverse  their  flow  direction  and  convect  out  of  the  computational  domain  along 
with  the  gas  flow.  If  the  initial  injection  velocity  was  greater  than  150cm/s,  we  observed  that  the 
droplets  were  able  to  change  the  gas  flow  field  in  the  plume  region.  This  disturbance  of  the  gases 
in  the  plume  region  resulted  in  generation  of  acoustic  waves  that  affected  the  flame  structure.  We 
note  that  the  droplet  even  for  these  high  injection  velocities  were  unable  to  reach  the  diffusion  the 
flame.  We  conclude  that  it  is  difficult  for  the  “Class  1”  droplet  sprays  to  penetrate  through  the 
plume  region  and  reach  the  diffusion  flame. 

Droplets  that  were  injected  away  from  the  plume  region  were  also  unable  to  overcome  the 
drag  force  exerted  by  the  cooler  co-flowing  air.  We  note  that  these  droplets  travel  further  down 
towards  the  diffusion  flame  than  the  droplet  that  were  injected  in  the  plume  region.  Since  the  top 
injection  boundary  is  approximately  30cm  above  the  burner  surface,  the  droplets  were  unable  to 
come  close  to  the  diffusion  flame  for  the  injection  velocities  that  we  have  tested.  Simulations  are 
being  performed  in  which  the  co-flow  air  velocity  is  set  to  zero.  Under  the  condition  of  zero  co-flow 
air  velocity  the  droplet  are  able  to  travel  downward  under  the  force  of  gravity  and  are  able  to  reach 
the  burner  surface.  The  droplets  were  found  to  entrain  into  the  diffusion  flame  resulting  in  flame 
suppression. 

4.  CONCLUSIONS 

Efficient  design  of  water-mist  systems  aims  at  obtaining  the  maximum  amount  of  sup¬ 
pression  for  the  minimum  amount  of  water-mist  added  to  the  system.  Numerical  simulations  have 
been  described  for  optimizing  the  various  mist  injection  characteristics  for  maximum  flame  sup¬ 
pression.  The  role  of  droplet  diameter,  injection  angle  (throw  angle),  mist  density  and  velocity  on 
entrainment  and  overall  flame  suppression  is  investigated  using  a  two-continuum,  time  dependent 
Eulerian  model.  The  net  suppression  ( AHj/AH0 )  for  various  water-mist  injection  configurations 
(discussed  in  the  Results  and  Discussion  section)  has  been  shown  in  Figure  14  as  a  function  of  mist 
to  fuel  flow  rate  ration.  Both  symmetric  and  asymmetric  spray  pattern  geometries  are  considered. 
This  figure  clearly  illustrates  that  net  flame  suppression  is  a  function  of  various  water-mist  injection 
characteristics.  The  goal  of  this  study  has  been  to  optimize  these  injection  characteristics  to  obtain 
the  maximum  flame  suppression  for  the  minimum  water-mist  to  fuel  flow  rate  ratio.  The  main 
conclusions  derived  from  this  study  can  be  summarized  as  follows: 

•  Under  base  injection  configuration,  the  water-mist  to  fuel  flow  rate  ratio  required  for  extinguish¬ 
ment  reduces  with  decreasing  droplet  size.  The  ratio  of  mist  to  fuel  flow  rate  (mmist/rhfuei) 

required  for  extinction  is  approximately  40  for  150/tm  droplets  and  reduces  to  10  for  50/m 
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diameter  droplets.  Flame  suppression  increases  with  higher  injection  density  or  lower  mist 
velocity  for  a  given  flow  rate  and  droplet  diameter. 

•  Under  base  injection  configuration,  net  suppression  reduces  as  the  spray  injection  velocity 
increases  (or  spray  injection  density  decreases)  for  a  given  mist  to  fuel  flow  rate  ratio.  Smaller 
droplet  diameters  produces  maximum  suppression  for  the  minimum  spray  injection  density. 

•  Under  symmetric  injection  (base  injection  of  water-mist  directed  at  45°  angle)  smaller  diameter 
droplets  produce  more  flame  suppression.  In  contrast,  under  asymmetric  injection,  larger 
diameter  droplets  produce  more  flame  suppression. 

•  Asymmetric  mist  injection  can  significantly  alter  the  flame  structure  and  flow  field.  In  gen¬ 
eral  asymmetric  conditions  produce  lower  suppression  even  with  higher  injection  densities  as 
compared  to  the  symmetric  case. 

•  Mist  injection  velocities  close  to  the  inlet  gas  velocity  result  in  net  suppression  that  is  nearly 
independent  of  injection  angle. 

•  Asymmetric  mist  injection  show  an  optimum  suppression  angle  of  45°. 

•  In  contrast  with  base  injection,  side  injection  configuration  clearly  shows  that  larger  diameter 
droplets  produce  more  flame  suppression.  Under  side  injection  configuration  flame  suppression 
increases  with  smaller  mist  density  or  higher  injection  velocity  for  a  given  flow  rate  and  droplet 
diameter. 

The  generality  of  these  conclusions  for  other  flames/fires  is  currently  being  investigated.  The 
following  additional  studies  are  being  conducted  to  obtain  a  detailed  understanding  of  the  physical 
process  involved  during  the  interaction  of  water-mist  and  flames.  These  results  will  be  presented 
in  subsequent  reports. 

•  Numerical  models  are  being  developed  to  simulate  the  burning  rates  of  two-dimensional  or  axi- 
symmetric  liquid  (methanol  or  heptane)  pool  fires.  An  energy  equation  is  solved  in  the  liquid 
pool  and  is  coupled  with  the  gas  phase  equations  through  appropriate  interphase  conditions. 
A  Clausius-Clapeyron  relationship  is  used  to  obtain  the  burning  rates  of  the  liquid  pool. 

•  The  temperature  profiles  above  a  liquid  pool  are  being  compared  with  experimental  data  ob¬ 
tained  over  a  similar  geometry.  The  burning  rates  of  liquid  pool  fires  is  also  being  compared 
with  experimental  data  available  in  the  literature. 

•  The  effect  of  the  various  water-mist  injection  characteristics  (base  injection,  side  injection  or 
top  injection)  on  overall  flame  suppression  is  being  investigated  for  methanol  pool  fires. 

•  The  changes  in  pool  burning  rates  as  a  function  of  the  various  water-mist  injection  character¬ 
istics  is  being  studied. 
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Fig.  1.  Schematic  diagram  of  a  2-Dimensional  diffusion  flame  burner  and  the  associated 
computational  domain  in  which  the  solutions  are  desired.  The  figure  also  shows  the  various  water- 
mist  injection  configurations  that  are  studied  in  this  report. 
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Fig.  2.  The  effect  of  droplet  injection  diameter  and  injection  velocity  on  overall  flame  suppression 
is  shown  as  a  function  of  mist  to  fuel  flow  rate.  The  arrows  in  the  insert  indicate  the  direction  in 
which  water-mist  is  injected. 


Fig.  3.  The  effect  of  droplet  injection  density  and  spray  pattern  (symmetric  or  asymmetric  injection) 
on  overall  flame  suppression  shown  as  a  function  of  mist  to  fuel  flow  rate.  The  arrows  in  the  inserts 
indicate  the  direction  in  which  water-mist  is  injected. 
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Fig.  4.  The  effect  of  varying  the  mist  injection  velocity  on  the  temperature  distribution  for 
asymmetric  mist  injection  at  45°  angle  to  the  air  co-flow.  The  initial  droplet  diamter  is  50^ 
and  the  injection  density  is  1000 drops/ cm3.  The  mist  injection  velocity  is  varied  from  0.25 m/s  to 
2.25m/s. 
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Fig.  6.  Temperature  contours  for  the  asymmetric  injection  configuration  of  water-mist  at  45°  angle 
to  the  air  co-flow.  The  water-mist  droplet  diameter  is  50 fi  and  the  injection  density  is  4000 
drops/ cm3.  The  mist  injection  velocity  is  varied  from  1.75ro/s  to  7.00 m/s. 


18 


Prasad  et  al 


Fig.  7.  Sectional  density  contours  above  a  methane-air  diffusion  flame  burner  for  asymmetrical 
injection  of  water-mist  at  45°  angle.  The  water-mist  droplet  diameter  is  50//,  injection  density  is 
4000  drops /cm3  and  mist  injection  velocity  is  1.75m/s. 
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Fig.  8.  Net  suppression  effect  of  water-mist  as  a  function  of  injection  angle.  Both  symmetric  and 
asymmetric  base  injection  configurations  are  considered. 
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Fig.  9.  Temperature  contours  during  symmetric  injection  configuration  of  water-mist  at  various 
injection  angles.  Droplet  diameter  is  150/i,  injection  density  =  100  drop/ err?  and  injection  velocity 
=  1.00m/ s. 
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Fig.  10.  Heat  release  rate  contours  (J/m3/s)  during  side  injection  of  water-mist.  Figure  shows  the 
exothermic  contours  (red)  at  the  flame  sheet  due  to  oxidation  of  the  fuel  molecule,  and  endothermic 
contours  (blue)  due  to  evaporation  of  water  droplets. 
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Fig.  11.  Velocity  vectors  {m/s)  superimposed  on  temperature  (K)  distributions  during  side  injection 

of  water-mist  into  a  methane-air  diffusion  flame. 
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Fig.  12.  Sectional  density  contours  above  a  methane-air  diffusion  flame  burner  for  side  injection  of 
water-mist.  The  water-mist  droplet  diameter  is  50/f,  injection  density  is  2000  drops/ cm3  and  mist 
injection  velocity  is  12.5m/s. 
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Fig.  13.  Net  suppression  effect  of  water-mist  as  a  function  of  mist  to  fuel  flow  rate  ratio  for  side 
injection  configuration.  The  effect  of  droplet  diameter  on  overall  flame  suppression  has  been  shown 
for  (asymmetric)  diffusion  flames.  The  arrows  in  the  insert  indicate  the  direction  in  which  water- 
mist  is  injected. 
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Fig.  14.  Net  suppression  effect  of  water-mist  as  a  function  of  mist  to  fuel  flow  rate  ratio  for  various 
injection  configuration  and  spray  pattern  geometries. 


